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The one dimensional Burgers equation in the inviscid limit with white noise initial condition is 
revisited. The one- and two-point distributions of the Burgers field as well as the related distributions 
of shocks are obtained in closed analytical forms. In particular, the large distance behavior of spatial 
correlations of the field is determined. Since higher order distributions factorize in terms of the one 
and two points functions, our analysis provides an explicit and complete statistical description of 
this problem. 



I. INTRODUCTION 

The Burgers equation for the velocity field u{x,t) 

d d d 2 

g^u(x, t) + u(x, t)—u(x, t) = v-^u{x, t) (1) 

has recently raised much interest because of its multiple connections to a variety of physical and mathematical 
problems. Background and references can be found for instance in the recent book |Q. The original Burgers problem 
concerned the statistics of the velocity field u(x, t) and of shock- waves in the inviscid limit v — > when the 
distribution of the initial velocity field u(x,0) is a (5-correlated Gaussian (white noise). It provides an over-simplified, 
but analytically tractable model of turbulence which has attracted a lot of studies over the last decades. In 
a considerable amount of work is done to calculate various moments of these distributions, but the distributions 
themselves were not obtained in closed forms due to the complexity of the analysis. The question has been addressed 
again by Tatsumi and Kida ||] and Kida Q. In ||, kinetic equations for the dynamics of shocks are used to derive 
scaling properties of the distributions, and in the second part of Q, Kida presents the result of numerical simulations 
for the distribution of the strength and the velocity of shocks. Recently, Avellaneda and E || and Avellaneda j(| 
have derived rigorous upper and lower bounds of the cubic type exp(— Cu 3 ) for the tails of the one-point distribution. 
Such cubic bounds have also been obtained in for the distribution of mass in the closely related problem of ballistic 
aggregation. 

In this work we revisit Burgers problem by providing closed analytical forms of the statistical distributions for the 
field and the associated distributions of shock- waves. Our main contributions are a simple formulae expressing the 
one-point distribution as integrals over the analytic continuation of the Airy function on the imaginary axis (formulae 
( |66| ) and (|7^)) as well as a detailed study of the clustering behavior of the two-point distributions. Since it is known 
that u(x, ijas a function of x (for fixed i) is a Markov process J^], the higher order distributions factorize in products 
of one and two-point functions. Hence our results give a complete solution to the one dimensional Burgers problem 
with initial white noise distributed data in the inviscid limit. 

In section II, we recall well known facts about the Burgers equation in the inviscid limit with the purpose to 
introduce the notations and the definitions of the one- and two-point distribution functions. In section III, using the 
notion of first hitting time, these distributions are expressed in terms of the basic propagator for a Brownian motion 
constrained by parabolic barriers. It appears that all statistical properties of the Burgers problem are embodied in 
the knowledge of three functions called here /, J and H . The functions / and J are calculated in section IV and the 
one-point distribution of fields and shocks are discussed. In particular an explicit formula for the strength distribution 
is obtained. These results have already been announced in [g| in the equivalent language of ballistic aggregation. The 
section V is entirely devoted to the study of the large distance behavior of the two-point distribution (the function 
H). Since the analysis is somewhat heavy, technical parts have been relegated in appendices. Finally the factorization 
of higher order distributions as well as their time dependence are discussed in the conclusion. 
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The situation considered by Burgers is particularly relevant to the non equilibrium statistical model of ballistic 
aggregation: it is known HQ that the dynamics of shocks in Burgers turbulence is closely related to the dynamics 
of the aggregating particles. White noise initial distribution of the Burgers velocity field corresponds to uncorrelated 
Maxwellian initial velocity distribution of the particles undergoing aggregation. Hence our results also solve this 
statistical mechanical model. A precise connection between the two problems can only be made in a proper scaling 
limit since ballistic aggregation always retains the discrete nature of particles whereas Burgers velocity field describes 
a continuous medium. This connection will be discussed in an other paper j^] . Notice also that the Burgers equation 
is equivalent to the KPZ equation |ic| l of surface growth processes with the height h of the surface given by u{x, t) = 
d x h(x,t). Burgers turbulence arising from other classes of stochastic initial data (see e.g. [pT[) or from the action of 
external random forces (see e.g. jl^l ), which is the subject of considerable current investigations, is not discussed in 
this paper. 



with 



II. GENERAL SETTING 

For convenience, we shortly recall the construction of solutions of the Burgers equation in the inviscid limit, see 
and references therein. Introducing the potential d^(x, t)/dx — u(x, t) together with the Hopf-Cole transformation 

V(x,t) = -2vm6{x,t), (2) 

one finds that the function 8(x,t) satisfies the linear diffusion equation 

d d 2 

It can be readily solved leading to the explicit solution 

u [ Xi t) = J ^ oc QO 1 v 2t/ K _ ;; (4) 

J^ rfyexp(-^F(a;,y,t)) 

where 

F(x,y,t) = ^^-^(y) (5) 

= -¥(», 0) = - f dy'u(y',0) (6) 
Jo 

which depends upon the initial condition. Burgers turbulence corresponds here to the situation where the initial 
velocity field u(x, 0) is a white noise process in space, or equivalently ip{y) is a two-sided Brownian motion pinned at 

V>(o) = o. 

In the inviscid limit v — > 0, the only contributions of the integrals in Eq.(Q) come from the minima of the function 
F(x,y,t), which depend on the initial condition through tp(y): 

£(x,t) = mm F(x,y,t) (7) 
y 

and we obtain 

u(x,t) = X ~^ t) . (8) 

Due to the scaling properties of the solution u(x,t), one can trivially take into account the time dependence of the 
problem. Indeed the scaled Brownian motion t a / 2 il>(yt~ a ) is equivalent in probability to i>(y) so that by (^) and (||) 
with a = 2/3 one has that £ 2 / 3 £ (7^73-, l) is equivalent to £(x, t) and t~ 1 / 3 u (jt/j, l) is equivalent to u(x, t). We study 
from now on the fixed time t = 1 solution u(x, 1) = u{x). It will then always be possible to recover the time-dependent 
solution through this scaling property as we shall sec in the concluding section. 

The minimum £(x) = £(2, 1) as a function of x can be found with the help of a nice geometrical interpretation 
of the solution. One considers a realization of the Brownian motion ip(y) and a parabola centered at x of equation 
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FIG. 1. Geometrical interpretation of the solution u(x) = x — £(x) for a given realization of the Brownian motion tp(y) 
which stays below a parabola of equation s(y) = (y — Xi) 2 /2 + C but on two contact points VKCO = (f» — Xi) 2 /2 + C and 
?/>(£i+i) = (£i+i — Xi) 2 /2 + C. A shock is located at Xi with strength fa = — £j and wavelength Vi — a — £i while 
r\i = /if/2 - /um- 



(x — y) 2 /2 + C (see Fig.|T|) and adjusts the constant C in order for the parabola to touch without ever crossing 
it. The coordinate of the contact point is the minimum £(x) leading thus to u{x) = x — £(x). Then one glides the 
parabola on the graph of ip(y) by a continuous change of its center x and C until it touches it for x = x\ on two 
contact points ^ and Thus at x = Xi, the function -F(x, y, 1) has two minima leading to a discontinuity of it(x), 
called a shock, where lim e ^o u(xi — e) = x\ — and lim e ^o tt(x, + e) = x\ — Ci+i- To make w(x) singled valued at a 
shock, we define it to be continuous from the left setting u{xi) = Xj — £j. 

A shock is characterized (see Fig.|l|) by its location Xj and two parameters which can be taken asQ 

Mi = 6+1 -& "strength", = Xi - & "wavelength". (9) 

Instead of i/j it will also be convenient to use the parameter 77^ 

Mi Mi % / lr .N 

»?< = -ir ~ VM, Vi = — . (10) 

The 

quantities of interest to be computed are on one hand the joint distribution densities p n {x\, ui\ X2, U2] ■ ■ ■ ;x„,u„) 
for the Burgers velocity field to have values in-between u\ and u\ + dui, ... , u n and u n + du n at points xi, . . . , x n , 



At time t, the strength is usually defined as the discontinuity fa/t = — ft)/* °f u(x,t) at a shock. 
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when average is taken over the realizations of the initial condition u(x, 0). On the other hand we will also consider the 
joint distribution densities of shocks p n (%i, Ml) Vii x 2i ^2, ??2! • • ■ ; x n , [i n , rj n ). We shall obtain the joint distribution for 
the Burgers velocity field u(x) from that of the variable x — £(x). At time t = 1, these two sets of variables coincide 
and we identify both distributions. 

Consider first the one-point distribution density p\(x,u) where u(x) — x — t;(x) = u. Because of translation 
invariance, pi(x,u) — pi(0,u) = pi(u) and u(0) = — £(0) = u. Hence p\{u) is the measure of the set of all Brownian 
paths ip(y) with ip(Q) = that have their first contact (f.c.)^] with a parabola y 2 /2 + C at £(0) = — u. As one can set 
the origin of coordinates to be at this contact point, it is given by the measure 

pi(«) =E{1>(y) < s u {y),y € R; f.c. with s u (y) at (0, 0)} (11) 

of the set of paths that stay below the parabola 

y 2 

Su{y) = -^-uy (12) 

and have their first contact with it at -0(0) = 0. By first contact in ([ll|), we mean that the path is strictly below the 
parabola tp(y) < s u (y) for y < 0, is assigned to pass at tp(0) = and is then such as t/Kz/) ^ s u(y) for y > 0. The 
expectation E{- ■ •} refers to Brownian paths running in the infinite "time" interval — 00 < y < 00. 

Likewise, the two-point joint density distribution P2(0, iii;x, Uz) = P2{x, 111,112) is the measure of the set of paths 
with ip(0) = that have a first contact with a parabola y 2 /2 + G\ (centered at the origin) at £(0) = —U\ and a 
first contact with a second parabola (y — x) 2 jl + C2 (centered at x) at £(x) = x — 112- Once again, we fix the origin 
at the contact point with the first parabola. Thus P2(x,Ui,U2) is the measure of the set of paths which stay below 
both the parabolas s^(y) = s ui (y) centered at u\ and a second parabola s^ 2 \y) centered at x + u\ of equation 
(y — x — ui) 2 /2 + C, while the paths have a first contact point ^(01 = with s^(y) and a first contact point 
ijj(x + ui~ U2) — q with s' 2 ) (y), where x > 0, x + u\ — > 0, see Fig. g []. In terms of this parameter q the equation 
of the second parabola is 

s^(y)= iy - x ; u ^ -^ + q (13) 

Now, q is arbitrary except for the constraints that the first contact point with s^\y) must be below the second 
parabola, namely s( 2 '(0) > 0, and that the first contact point with s^(y) must be below the first parabola, namely 
+ u\ — U2) > q- This leads to the condition —q% < q < (72 with 

qi = q 1 (x,u 1 ,u 2 ) = -[x + U\ - u 2 )(x + Ui + u 2 ), q 2 = q2(x, Ui, u 2 ) = ^(x + Ui - u 2 )(x - Ui - u 2 ). (14) 
Hence 

P2 (x,u l7 u 2 ) = E {^(y) < s {1) {y), i/)(y) < s^iy), y G R; 

f.c. with s (1) (y) at (0,0), f.c. with s (2) (y) at (x + m- u 2 ,q), -qi < q < q 2 } • (15) 

The distributions p±(ui) and p 2 (x, Ui, u 2 ) have the normalizations 



and 



duxpx{ux) = l (16) 



du 2 p 2 (x, 111, u 2 ) — pi(ui), lim p 2 (x, m, u 2 ) — 6(ui - u 2 ) (17) 

x— >0 



2 Consideration of the first contact (or hitting) point is consistent with the left continuity of u(x). If there is a shock at Xi, 
u(xi) — linie^o u(xi + e) = fi+i — > 0, implying that £j has to be the first contact with the parabola. 
3 The case x + ui — 112 — 0, i.e. when the two contact points coincide, is discussed in the next section. 
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FIG. 2. Brownian interpretation of the two-point distribution of the velocity field P2(x, Mi, 112). The Brownian paths stays 
under the parabolas s^'ty) = y 2 /2 — u\y and s' 2 '(t/) = (y — x ~ u±) 2 /2 — wl/2 + q but on two contact points 4>(0) = and 
ij)(x + iii — U2) — q where —(as + ui — U2)(a; + ui + U2)/2 < q < (x + ui — U2)(x — ui — U2)/2. 

The distribution of shocks are defined in the same manner. By translation invariance p\ (x, /j, 77) = Pi(0,/i, 77) = pi(n, 77) 
is independent of x. It is given by the measure of the set of paths that have two contactsjj with the parabola 
s v (y) = y 2 /2 — vy (recall that v = ^ — a first contact at V'(O) = and a last contact (I.e.) at — 77, (see Fig. 

/Oi(M,?y) =E{i/;(y) < s„(y), y G R; f.c. with s„(y) at (0,0); I.e. with s v (y) at (//, 77)} . (18) 

The joint distribution p2(0, /Ui, 771; x, /i2, 772) of two shocks at distance a; is the set of paths that have two contacts 
with the parabola s Ul (y) as above and two contacts with an other parabola whose characteristics will be given in the 
next section. Notice that the centers of the two parabolas are separated by a distance x. 

All quantities will be eventually expressed in terms of the transition probability kernel for Brownian motion in 
presence of parabolic absorbing barriers |l3|Ji"4| ]. Consider the conditional probability density K v {n\, 771, (/,%, 772) for 
the Brownian motion tp(y), starting from ip(pi) — 771, to end at ^{^2) = 772 while staying under the barrier ip(y) < 
8v{y) = y 2 /2 - vy for m < y < fi 2 

^1/(^1,771,^2, m) = Em, vi Wv) < s ^(y)' Ml < V < ^2; i>(H2) = 772} • (19) 

It thus satisfies the diffusion equation 

d^Kvbn, 771,^2,772) = -^d^ 2 K u (fj,i,r]i, fi2,V2) (20) 

with K v (ii, 771 ,^,772) = 5(771-772) and K v ([ii, s u ([ii), /i 2 , 772) = K v (pi, 771, /x 2 , s„(/i 2 )) = 0. The parameter D in Eq.@) 
characterizes the initial condition. To solve this equation it is convenient to consider the shifted stochastic process 
4>{y) = if>(y) — s v {y) which is a Brownian motion with a parabolic drift. Clearly 

K u (fll, 77l,/X2,772) = K(jll,r]l - S v (fJ,l),fjL2,T)2 - S V (H2)) (21) 



4 If a path has more than two contacts with the parabola, the shock parameters are obtained in terms of the coordinates of 
the first and the last contacts. 
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where K satisfies the diffusion equation with drift 



D 



(22) 



with K(n,(f>i, /J,,4> 2 ) = S(4>i — (j)2) and Dirichlet boundary conditions K {jx\, 0, // 2 , 02) = K(fj,\, (f>\, fj, 2 , 0) = 0. Eq.(|2 
can be reduced to a diffusion equation with linear potential by the transformation 



G(/Zi,0i,/x 2 ,0 2 ) = -ftT(^i,0i,/i 2 ,0 2 )exp 



Then the propagator G is the solution of the equation 

d D d 2 1 



dfj, 2 2 d(f> 2 2 D 



02S"(^2) G(/il,0i,/Z 2 ,0 2 ) = 0, 01, 2 < 



(23) 



(24) 



d a 2 



with G(/i, 0i, (i, 4>2) = — 4>2) and Dirichlet boundary conditions at the origin, G(/^i, 0, /i 2 , 02 ) = G(/ii, 0i, /i 2 , 0) = 
0. Since s"(/z) = 1, this equation can be solved with the help of the spectral decomposition of the operator — y 
i0 2 leading to 



G(/il,01,/i2,02) = ( ) X] 



1/3 



fe>l 



Ai(-(2/^ 2 ) 1 / 3 1 - a; fc )Ai(-(2/ J 2 ) 1 / 3 2 - c^) 

(Ai'(--,)) 2 • ( ] 



The Airy function Ai(iu) p5), solution of 



/"H - wf(w) = 0, 



(26) 



is analytic in the complex w plane, and has an infinite countable numbers of zeros —LUk on the negative real axis, 

< U)\ < LU2 < ■ ■ ■■ 

Finally coming back to K v with the help of and introducing the explicit form (|12j) of s v (y) leads to 



K u (^!,m^2,V2) = G(^i,0(^l)),Ai 2 ,0(/i2)) 

1 



x exp 



0(Mi)(Mi 0(M2)(M2 -v) + 



(Pi - vf (A*2 - v) 



3 \ 1 



£> ' Mr " v " ' 6 6 

with 0(jui) = r\\ - 8 v ([ii), 0(/i 2 ) = m - s v (fi2). Note the symmetry K„(fii, 771, /x 2 , %) = #-*(-/«2, »?2, 



(27) 



III. DISTRIBUTIONS AND TRANSITION KERNEL 

In this section we relate the distribution functions to the transition kernel K u (fii, rji, // 2 , T) 2 )- We first treat the case 
of a single first contact point by computing the (conditional) probability density 



E flltVl {'>p(y)<s u (y),iJ, 1 <y<H2; f.c. with s^y) at (//, s „(//)) ; ^(// 2 ) = r? 2 } 



(28) 



that a Brownian motion starting at point (/ii,7yi) ends at (/U 2 ,i7 2 ) while staying below the parabola s v {y), and has a 
first contact point "0(^0 = s iAy) a t "time" /i, with /ii < fx < /i 2 and 771 < s„(/ii) and ?7 2 < s u (fj,2). This enables us to 
write the probability pi(u) (|ll|), where the expectation is taken on paths that run in the whole "time" interval y € R, 



as 



Pi(it) = lim lim / drjx I dr\2 

A11 ^-oo Al2 ^oo i /_ oc J_ ca 

E Hi,viWv) ^ s «(2/)ii"i < 2/ < A*2; f.c. withs u (y) at (0,0); t/K/^) = ^2} 
As in the preceding section it is convenient to consider 

E ^,4>A4>{v) < 0, /ii < y < /i 2 , f.c. with the origin at (^,0); 0(/j 2 ) = 02} = -d^P^^^^ifi) 



(29) 
(30) 
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the quantity corresponding to ( |28| ) for the shifted process 4>{y) = ip(y) ~ s ^(y)- It i s the (conditional) probability 
that a drifted Brownian motion </>(/i), starting at point (/ii,^>i), ends at ((J- 2 ,4> 2 ), stays negative </>(y) < and has a 
first contact with the origin at "time" \i (</>(/-*) = 0)- The desired quantity ( p8| ) is obtained by setting tfii — 0(/ii) = 
Vi ~ Sv{^i)t4>2 = ^(1^2) — V2 ~ 81/(^2) m (pdj). We have also written that (30) is the density of the probability 
(/i) that, under the same constraints, the path has its first contact with the origin at some "time" larger 



P 



Ml,?>i;M2,<?>2 



or equal to /1. This probability is given by (for basic notions on first hitting time see |lq| ) 

f° — — — 



(31) 



Indeed one considers the paths starting from (fj,i, <f>i) that stay negative up to (fi, <p) and then vanish at some "time" 
larger or equal to \x. The probability density for the later part is given by the measure of paths staying below the 
displaced barrier <j)(y) < e diminished by that of paths staying below the origin <fi(y) < as £ -> 0, namely by 



lim-(K(p, 

e^O e 



K(/J,,(f),lJ,2,<h)) 



(32) 



This leads to (|l|). Introducing ( |3"l| ) in (^0|) and using the forward diffusion equation (^2|) as well as its backward 
equivalent, we find after several integrations by parts that 



E iJ.x,<t>i{^(y) < °. Mi < V < M2; f-c with the origin at (/j,0); <f>(fi 2 ) = $2} 



D 



d^K(pi, 0i, /i, 4>)d^K(tJL, (f>, M2, 



(33) 



0=0 



Coming back to the original variables our probability (g8() reads 

r. 

-d v K v (pi ,r\i,\x, ri)d v K v ([i,r\,[L 2 , 772 



t^mi^iv) ^ s Av)^ Mx < V < M2; f-c withs^y) at (p,, s v (p)); ip(fx 2 ) = m) 

D 
~ ~2 



(34) 



with K v given by Eq.@. 

Let us define the function J(v) to be 



ID 



J( v ) = -\hr lim 

Z M2— >00 



drj^K ^(0,7], fi 2l r]2) 



(35) 



Then, from ( pi] ) and ( j29| ) it is now straightforward to find the expression of the one-point distribution of the velocity 
field 



Pi(u) = J(~u)J(u) 



(36) 



where we used the fact that -/^(/Ji, 771, 0, 77) = if_„(0, 77, — fii, 771). 

We come now to the two-point function ( |l5| ) which involves a first contact at y = with the parabola (y) and 
a first contact at 7/ = x + 7i! — u 2 with the second parabola s <y2 \y). We consider first the situation where these two 
contact points are distincts, i.e. when the strict inequality x + u\ — u 2 > holds. Since u(x) has slope equal to one 
except at the location of shocks, this corresponds to velocity fields with tj(0) = ui,u(x) = u 2 that have at least one 
shock in the interval [0, x). 

When x + ui — u 2 > 0, each contact gives rise to an expression of the form (|34|) with appropriate parameters 
(see Fig. ||). The first contact with the parabola s^(y) — s Ul (y) is as before. After the "time" p* (coordinate 
of the parabolas intersection s^(fj,*) = s( 2 )(/j*)), the paths are found under the second parabola s^ 2 \y) ( |l3| ) with 
corresponding propagator K s {2) (p\, r\\, \i 2 , rj 2 ) and first contact at (x + ui — u 2 , q). The corresponding probability is 
given by the following arrangement 



D 
~2 



d v K Ul (m , 771, 0, ii)d v K Ul (0, 77, [i* , 77') 



r/=Q 



D , 



XyS^is) (n*,r)',x + ui - u 2 ,T])d v K s (2)(x +ui- 7/2,77,^2,772) 



(37) 
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which has to be integrated on T)i, rj 2l if and q in the appropriate ranges and taken in the limits \i\ — > —00, /12 — > 00. 
The propagator associated with the second parabola can be written in a coordinate system where the second contact 
point is again located at the origin, namely 

K s (2)(ni, 771, ^2,m) = K U2 (ni - x- ui +U2,J?i - q,H2 ~ x - Ux + u 2 ,r] 2 - q) (38) 

One finds finally 

p 2 (u 1 ,U 2 ,x) = J(-Ui)H(x,Ui,U 2 )J(ll 2 ), x + u 1 -u 2 >0 (39) 

where the function H(x, v\ 1 v 2 ) is defined as 



H(x,v u v 2 ) = J d 1 I <V'VC, (0. '/■/''•'/) 



d v K V2 (/x* - x - V\ + v 2 , rj - q, 0, 77) 

77=0 



(40) 



7J = 



The integration limits gi = (71 (a;, z/i, v 2 ) and (72 = 92(2;, V\i v 2 ) are given by (|14]). The intersection point between the 
two parabolas has coordinate (/i*, 77*) = ((g + qx)/x, /i* 2 /2 — uifj,*). 

We now determine the contribution to ^2(^1 ^11^2) of the set of velocity fields u(x) that have no shocks in [Q, x) 
(i.e. when x + u\ — u 2 = 0) with the help of the normalization (p"7|). The set of Burgers fields with u(0) = Ui can 
be divided into the union of two disjoint sets, those having at least one shock in [0, x) and those having no shocks 
in [0, x). As seen before the first set corresponds to Brownian paths having two distinct contact points and from the 
previous discussion its measure is given by J" 1 31 du 2 J(—ux)H(x, u\, u 2 ) J(u 2 ). The second set corresponds to the 
case x + Ui — u 2 = when Brownian paths have a first contact point tp(0) — at the intersection of the two parabolas 
s Ul (y) and s Ul+x< (y) with measure 

E{ip(y) < s Ul (y), 77 < 0; t/%) < s Ul+x (y), y > 0; f.c. with s Ul (y) at (0, 0)} = J(-«i) J{u x + x) (41) 

The result ( pi] ) is derived by a slight extension of the calculation that led to ([36]) . The measures of these two sets 
sum up to pi{u\) 

/Ul+X 
du 2 J{~ui)H(x,ui,u 2 )J(u 2 ) —pi(ui). (42) 
-00 

Hence we conclude from ( [l7| ) that the complete form of p%{x, Ui, u 2 ) is 

p 2 (x, ui,u 2 ) = J(—Ui) [S(x + ui - u 2 ) + 9(x + 7/1 - u 2 )H(x, ui,u 2 )] J(u 2 ). (43) 

A quantity of interest is the probability density P[o }X )( u i) f° r the Burgers field to take the value u x at x = while 
there is no shock in the interval [0,x), i.e. u(x) = u\ + x. This is precisely the quantity (^lj), namely integrating ( |4^ ) 
on 7i2 with H omitted 

P[o,x)(ui) = J{-u x )J{ui +x) (44) 

and thus 



P[o,x) = / duiJ(-ui)J(ui+ x) (45) 
J — 00 

is the distribution of intervals of length x without shocks. 

We turn now to the shocks distribution functions. According to the discussion of previous section Eq. (p"8|), we use 
Eq. ( |34| ) to write the one-shock distribution function considered as a function of the parameters /1, 77 (see Fig. [j]) 



Pl ( M , 77) = J{-v)i{^ V )j(-n + v) = j { - f - + j j n) J ( -£ - -i ) 1 4(i ! 



where the function /(tx, 77) is defined as 

n 

(47) 
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The two-shocks distribution p 2 (0, /Zi, 771 ; x, M2j ^2) = /°2(^; Mi? 171 ; , ^2) (considered as a function of the shock 
parameters 771 , /ii and 772, M2) can be written as 

pi{x, fix, 771, ^2,172) = J(-^i)/(/xi,?7i) [<5(a? + r^i — i/ 2 — Mi) 

+6»(a; + z^i - z/ 2 - fix)H(x, -pi + ^1,^2)] I(P2,V2)J(-P2 + v 2 ) (48) 

with i/j = /ii/2 — rji/ fii, i = 1,2, and the functions /, J and i? are defined above. 

We denote p 2 nn \x; fix, 771; /i2, 772) the probability density of two nearest neighbors shocks separated by a distance 
x; pj""' ( x ! Mi) ^lj M2j ^2) is given by the formula ( |48| ) with the 7J function omitted. Then the conditional probability 
density p( nn >(fj,x, Tix\x, fi%, 772) that given a shock fix, T)x &t x = 0, the next shock P2, r] 2 occurs at x > is found to be 

(nn), I \ P2 ln) (x- 1 fli,7 1 i- 1 fl 2 ,V2) 

Mi(Mi>»7i) 

V Mi M2 2y 
This conditional probability has the normalization 

1 />oo />oo 

da;/ dfi 2 dif 2 p {nn) (fii, m\ x , M2, m) = 1 (50) 



/• OC 

io ' 



'0 

which leads to the following integral relation between the functions / and J 

m - r * /> • (t - 1 - ") /( "- " >j B-5) <5i) 

This analysis shows that the one-point and the two-point distribution functions of the Burgers velocity field u(x) 
as well as of the statistics of shocks are entirely determined by the knowledge of three functions /, J and H defined 
in (47), (P^) and (|40|). Finally, these last three functions can be computed from the basic transition kernel K u given 
by Eq.(gT 



IV. THE FUNCTIONS I AND J AND THE ONE-POINT DISTRIBUTION 



In this section we giv e explicit expressions for the functions I(fi,r)) and J{v) defined respectively by Eqs.([47|) and 
([35]). Through Eqs. (1361 J46| ) , we will then obtain explicit forms for the one-point distribution function of the velocity 
field px(u) and of the shocks 77). 

Using the form (p7|) of the transition density K v in Eq.(|47|) we have that 



I(fi, 77) = 2a 3 exp 



!L 
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j(m). 



We set a = (2ZTT 1 / 3 and 



(52) 



(53) 



k>l 



where — u>k, k > 1, are the zeroes of the Airy function. This last expression has already been found by Burgers JJ2[. 
Our point here is to give a closed form for the function J{v) and thus for the one-point distributions. Inserting ( |27| ) 
in (35) and changing the variable x = — (2/D 2 ) 1 / 3 (rj2 — s !/ (m2)) leads to 



J(iz) 



r a hm e~ a 3 [(a-- ) 3 +- 3 l/3 

/i — t-OO 



/•DO 

1/0 k>i 



Ai(x-w fe ) 



Ai'(-wfc) ' 

It is convenient to introduce the following integral representation of the sum for p > 

, Ai(w + x) 



fc>i 



A\{x-u k ) 1 j d ^ caw[l Ai(w + x) 
Ai'(-u k ) 2m Jc Ai(w) 



1 

2nd 



dwe aw >*- 



Ai(w) 



(54) 



(55) 
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where the contour C runs just above and below the negative real w— axis encircling the zeros of the Airy function. 
From the asymptotics 

Ai(io) = (4iry/^y 1/2 e- 2w3/2/3 (l + C(w" 3/2 )) , \w\ -> oo, |argw| < ir (56) 

one deduces that for w — \w\e lS , n/2 < 6 < n, one has |Ai(iu + xMAi(w)\ ~ exp (— xlwl 1 / 2 cos(0/2)), |to| — > oo, 
cos(#/2) > 0. For 9 ~ tt, the factor e ™^ ensures the convergence in (p5j). Hence for \i > one can deform the contour 
C and show that the unique contribution to the integral comes from the imaginary axis —ioo < w < ioo leading to 
the last part of the identity (p5|). After exchange of the integrations order one finds 

-i />zoo „a[iw /* oo 

J(v) = y/a~ lim e" a [ ^- v) +v ]/3 / dw—— / dxe- aa:(l/ -^Ai(x + «;). (57) 

A»-oo 2l7T 7_ JOO Al(l0) 7q 

To proceed we determine first the Laplace transform of f(x) = Ai(x + w), w fixed, 

/>oo 

f(s)= / dxe- xs f{x). (58) 



JO 

The function /(x) is solution of the second-order differential equation 

f"(x)-(x + w)f(x) = (59) 
with /(0) = Ai(w) and /'(0) = Ai'(w). The Laplace transform of this equation is 

/'(*) + (s 2 - uj)f(s) = sf(0) + /'(0). (60) 



with solution 



M = (f(0) + J' da ((7/(0) + /'(0))c— +^/ 3 ) e ws ~ s3 / 3 (61) 



/(0) = / dx Ai(x + w) = -tt [Ai'{w)Gi(w) - Ai(w)Gi'(w)] (62) 



where Gi(iu) = tt" 1 J™ dt sin(i 3 /3 + wt) (Tl. 

Inserting this Laplace transform in Eq (|57|) and using various properties of the Airy functions [l5| leading to the 
identity 

/(0) - f da (af(0) + /'(0))e-" CT +- 3 / 3 = 1, (63) 



we eventually find 



with 



J(u) = Va~e~ a3 » 3 / 3 J(is) (64) 



= — / dw——. (65) 
2trr J_ ioo Ai(w) 

Note that this integral is convergent for positive and negative v. 

With the explicit form (|64| ) of the function J{v), the one-point distribution function p\(u) of the velocity field is 
given by 

Pl (u) = J{u)J{-u) = aJ{u)J{-u) (66) 

which is plotted on Fig. ||. 

Defining the moments of the distribution as (it") = J duu n pi(u) we have (u) = as Pi(u) — Pi(—u) and (u 2 ) = 
toi(Z)/2) 2 / 3 with a constant mi ~ 1.054. The normalization ( |l6| ) is verified as, from (|66|), 
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FIG. 3. The one-point distribution function (p6|) for the velocity field pi(u) as a function of u for D = 1/2 (a = 1). 



r°° i r 

/ dup!(u) = — 

J-oo 2l7T J_, 



dw 



[Ai(wW 



(67) 



can be shown to be equal to one. 

To determine the asymptotic behaviour of pi(u), we remark that for positive u, we can close the contour in (|65| ) to 
encircle the poles of the integrand and thus express S(u) as a sum on the zeros of the Airy function 



u>0 



(68) 



Hence J(u) ~ e allWl /Ai'(— Wi) as u — > oo. The behavior of J(u) for u — > — oo can be determined with Laplace 
method to be l7(u) ~ — 2au exp(a 3 u 3 /3) and so the large |u| behavior of pi(u) reads 



2a 2 M 
PlW ~ TT77 7 ex P 

Ai (-wi) 



a 3 |w| 3 



a|u|a;i 



(69) 



This result is of course compatible with the bounds found in Th. 1 of Q, but cubic bounds cannot be saturated 
because of the additional exponential decay exp(— a|u|o;i). It is interesting to remark that, starting form a Gaussian 
distributed initial velocity field u(x, 0), the field immediately evolves to a distribution which is not Gaussian but 
behaves as Eq. (po]) . 

Let us turn now to the one-shock distribution function pi(fi,rj). Collecting results Eqs. (E^,B2lp3), we find 



fi 2 



fi 2 



(70) 



with X and J defined in (|53|) and (|65|) , respectively. 

One can compute the shock strength distribution defined as 



Inserting 



Pi(p) = / drtpi{fi,Tf). 
in this last equation, we find after the change of variables w = iQ and rj — at]/ /i 



(71) 
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FIG. 4. Shock strength distribution for D = 1/2, (a = 1 in Eq.(ph 



= 2a 3 /iX(it) 



dCi 



which reduces to 



with 



(2^)2 y^^V-oo^Ai^COAi^) 



Pl (fi) = 2a 3 fiX{n)n(fj,) 



-ia M (Ci+C 2 )/2 /-oo 



1 

2iTT 



dw- 



(72) 
(73) 
(74) 



The form of the shock strength distribution (|73j) is plotted on Fig. W Notice that L (/z) is the space covered by the 
shock strength in a box of size L; it is equal to L and one has thus d/i /zpi(/i) = 1. 

We can now determine the behavior of the shock strength distribution for small and large shocks. For 0</i<l, 
we use the normalization condition ( |l6| ) to find TL(n) = 1 + 0{p) while the behavior of X(/i) can be determined 
from the large k asymptotic behavior of the zeros of the Airy function ujk — (37rfc/2) 2 / 3 + 0(k^ 1 ^ 3 ) to give X(/x) ~ 
(2 v ^ : (a/Li) 3/2 )" 1 . One thus get 



Pi(m) 



+ G(At 1/2 ), M^O. 



(75) 



The divergence /i -1 / 2 , as /i — * 0, has been found in |5) and seen in numerical simulations 

On the other hand, for large /i, one can estimate the behavior of the function 7i(/i) by the Laplace method to find 
7i(/i) ~ (ira 3 /! 3 ) 1 / 2 exp(— a 3 /i 3 /12). The behavior of the function X(ju) is immediately given by the largest zero of the 
Airy function to give I(/i) ~ exp(— u>id/z). We thus have 



/ fl 3 3 \ 

pi (it) = 2 VW /2 /i 5/2 exp f jj- - uiafi J , 



(76) 



Let us consider now the shocks wavelength distribution. The one shock distribution (|70|) can be written for the 
strength-wavelength variables (/i, v) as^] 



5 The additional fj, factor is the Jacobian of the transformation (n,r/) to (/x, v). 
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FIG. 5. Shock wavelength distribution pi(v) for D = 1/2, (a = 1 in Eq.(||)). 



pi([i, v) = 2a i nJ(-v)l(n)J(v - fj,). 



(77) 



Considering the variable v' — u—fi/2 we find that the (fi, v') distribution is symmetric in v 1 , implying (y 1 ) = and thus 
(v) = = |. The wavelength distribution pi(z/) = / °° dfj, pi(fi, v) is plotted on Fig. |[ Its asymptotic behaviour is 
found to be p\{y) ~ C+v 3 exp (— a 3 v 3 /3 — avuoi) , f — * oo, and p\{v) ~ C_ exp (a 3 ^ 3 /3 + ai^i) , v — * — oo. Remark 
that the wavelength distribution is not symmetrical around v = i. 



The density distribution P[o yX ) of intervals of size x with no shocks (45) is given by 

/oo 
(foil J(-ui)J(x + lii) 
-OO 

1 



7T 

ox 6XP V 12 J (2tt^ 2 



joo exp 
duj2 

■ioo 



Ai(wi)Ai(w 2 ) 



(78) 



which is plotted on Fig. [| Since lim :E _ > oP[o,2;)( ,i i) — Pi( u i)> ( see Eq.(^4])), and p\{u) is normalized (|l6|), we have 
lim 2 .^o?'[o.x) — 1- Asymptotically we have for x — > 00 



P[0,x) 



— exp(-4f 



aujix 



ax [Ai'(-wi)] 4 



l + 0[ - 

x 



(79) 



V. CORRELATIONS 

In this section we study the two-point distributions of the Burgers velocity field and of the shocks in the asymptotic 
limit x — > 00, keeping all the other arguments fixed. From (B3J) and (Eq), we have for x large enough 



and 



P2(x,ui,u 2 ) = J{-ui)H(x,ui,u 2 )J(u 2 ) 
P2(x,/j,i,r)i, M2,%) = J{-V\)I{p-i,'ni)H{x,-pLi +Vl,lS2)I(p2,ri2)J(-H2 + V2) 



(80) 
(81) 
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FIG. 6. Distribution P[o,x) of intervals [0,x) which contains no shocks for D = 1/2 (a = 1 in Eq.(|78[)) 



with the function J and / given by Eqs. ( p4| , p2|) and where the function H is defined by Eq.(|4C]) with Vi = (ii/2 — r\ij /ii 
Our main results are 



p 2 (ux,U2,x) - px(ux)pi(u 2 ) 
8^ 



■ exp 



a 3 x 3 



a i/2 x 5/2-*-y | 2 
and similarly for the distribution of shocks 



aujix ) exp (— auji(ui — u 2 )) J{—ui)J(u 2 ), x 



(82) 



32 /7F 

P2(mi, m,V2,m, x ) - Px{^Wni)pi{^x,r}{) a 11/2 ^ exp 



auj\x 



x 5/2 -r y 12 

x exp (—awi (i/x - ^2 - Ml)) J {- v i)1{ni)1{H2)J {- P2 + v 2 ), x — » 00. 



(83) 



Wee see that long distance correlations are very weak since they are again dominated by the cubic decaying factor 
exp(— a 3 a; 3 /12). 

Clearly, in view of © and ©, this asymptotic behavior is determined by that of the function H(x, v%, v 2 ). First 
we write H(x 1 vi 1 v 2 ) in explicit form by introducing J27| ) in (40). It is useful to remember that by the definition of 
(p*,rj*) one has if = s yi (/i*) = q + Sj/ 2 (/i* — x — v\ + v 2 ). To bring the expression in the most symmetric form the 
change of integration variables 



C = 0D 2 /2) 1/3 (rf -rf), 
1 

Jx 



1 2 



< C < 00 

X X 

■y/xri <r< y/xr 2 , n = - + v\, r 2 — - - v 2 



(84) 



turns out to be adequate. Then, with a = (2D) x / 3 , one obtains 
H(x,vx,v 2 ) — 2a 3 exp 



, —arv? + a 3 is 9 , ,- 

■" ! 1 \ ,|-CX|) 



ax 



d(e a<:x ex P 



fci ,fc 2 



-aaj kl [ n + —= 



12 

auj k2 r 2 



s/xr 2 



dr exp (— a 3 r 2 ) 
Ai'(-c fel )Ai'(-a; fc2 ) ' 



(85) 



Our main concern is to determine the asymptotic behavior of this expression as x — > 00. We give here the main steps 
of the calculation while details and justifications are given in the appendices. 
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To get the basic clustering properties of the model, we expect that lim^oo H(x, f±, v 2 ) — J{v{) J{— v-i) with J(y) 
given by the integral in the complex plane Eq.(0). It is therefore natural to replace the sums on the zeros of Airy 
functions in j85| ) by appropriate contour integrals, as in Sec. IV, 

Ai(C~w fel )Ai(C-^ fc2 ) _ 



ex p 

fel ,&2 

1 



-au kl \ri + -fi 



/ dw\ I dw 2 exp 








awi fri 4 




+aw2 ( r2 ~^)_ 


lc x JCx 









Ai(C + wi)Ai(C + w 2 ) 

Pd{wi)k\(w2) 



(86) 



For a given x the contour C x is chosen as the parabola with branches w (p) = —p ± iax^/p, < p < oo. This 



awi [ n + ^ ) + aw 2 \r 2 ~ — 



, 5Rwi < 0, dtw 2 < 



contour will be convenient to determine the large x asymptotics of H(x,i>i,i> 2 ). The integrals (86) on C x converge 
for r fixed because of the exponentially decreasing factors exp 

°> Tl + 72 > °> r2 ~ 75 > ( see A PP endix 0)- 

Next we exchange the C-integral with the contour integrals to obtain 

\fx exp 



H(x, v\, v 2 ) = 2a exp 
1 



12 



dr exp i—a r 



3^,2^1 



(27Ti) 2 



where 



dwi / du> 2 exp 



B(x,wi,w 2 ) 



awi r 1 + 



aw 2 r 2 



B(ax, Wi, w 2 ) 
Ai(u>i)Ai(w 2 ) 



dCe Cx Ai(C + wi)Ai(C + w 2 ) 



(87) 



(88) 



is the Laplace transform of a product of Airy functions evaluated at the negative argument —x. This Laplace transform 
is comp uted in Appendix H and is given as the difference of two terms B(x,wi,w 2 ) — Bi(x,wi,w 2 ) — B 2 (x,wi,w 2 ) 
(see Eq.(pS|)). We set H(x) = H 1 (x) — H 2 (x) with H 1 (x) (resp. H 2 (x)) the contribution to ( ]87| ) of B 1 (x,Wi,w 2 ) 
(resp. B 2 (x,Wi,w 2 )). Then 



H^x) 



,5/2 



■ exp 



dr 



1 



(2«) 2 



dwi / dw 2 h 1 (r/,w 1 ,w 2 ) 



with 



exp 



wi, w 2 ) 



- 3 --w; 



V 



+ aw\V\ — aw 2 v 2 



(89) 



(90) 



Ai(wi)Ai(w 2 ) 

It is shown in appendix ^| that the multiple integral in (|89| ) is absolutely convergent. As x — > oo, the contour C 2 - 
eventually opens to the imaginary axis of the w-plane. Hence one sees (formally) on (|8^) that 

exp (aw\V\ — aw 2 v 2 ) 



lim Hi(x) — a exp 



-a 3 j/ 3 + a 3 ^| 



dwi 



dw-2 



Ai(wi)Ai(u> 2 ) 



= J{yx)J{-Vi) 



(91) 



where the function J(y) is defined by Eq.(J6J). More precisely one finds that the asymptotic behavior of H\(x) is 
given by (Appendix O 



Hi(x)= J(v!)J(-v 2 )+0 exp 



12 



-(1 + c) 



c> 0. 



(92) 



Inserting the expression B 2 (x) (B9) in (|87| ) one finds 



ff 2 (x) = 2a 5/2 exp 
1 



-a 3 ^ 3 + a 3 ^| 



dr e 



(27Ti) 2 

x exp 



dwi / dw 2 exp awi ri + — + aw 2 r 2 — / dy^fy 



y 

12 



101 + W 2 (l£?l — ?« 2 ) 2 

[y - ax) 



1 1 

ax y 



g(y,w 1 ,w 2 ) 
Ai(uii)Ai(w 2 ) 



(93) 
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Because of the convergence factors exp (^awi (n + ~^j= \ + aw 2 (r 2 — ^ e con t° urs Cx can be closed and the 

corresponding integrals can again be eval uated at the zeros of the Airy functions (the arguments are similar to those 
given in appendix 0). Then the relation (B1C) permits to simplify the result to 



H 2 (x) = 2a 5/2 exp [ L — ) / dre " 



— y/xrx 



X ^2 ex P I [ r J 



auj k2 r 2 



x 



2 4 \ax y t 

To compute the large x behavior it is convenient to make the change of integration variable y = + ax giving 

'("*) 



(94) 



H 2 (x) = 2a exp 



-a 3 v 3 + a 3 v 3 \ ex P ( 



c 5/2 



-G(x) 



(95) 



with 



G{x) = dz- 
Jo 



exp 



\ Tlx® 



az a z 

4x a 4 



\fxr 2 



dre ° 3r ^2 ex P ( _ ( w fci _ ^fc 2 ) 



Aa(zx + ax 4 ) 



•y/x 2air 2 

Letting formally a; — * oo on this formula gives the asymptotic behavior (details are found in appendix [d]) 



8V^F (-a 3 vl + a 3 vl 
H ^ x ) ~ a 3/ 2a .5/2 ex P ( 3 tuji(i/i - i*) J exp 



aui\x 



(96) 



(97) 



where — u>\ is the first zero of the Airy function. 

Inserting the asymptotics Eqs. (|97| , |92] ) in the expression for the two-point distributions Eqs.(^,|l|) leads to the 
results Eqs.tfll!). 



VI. CONCLUSION 



To conclude, we remark that the previous results allow for a complete statistical description of the Burgers field. As 
mentioned in the introduction, for white noise initial data, u(x) is a Markov process as function of x Thus with 
P(x 2 ,u 2 \xi,ui) — p 2 (x 2 — xi,ui,u 2 )/pi(u\) the transition kernel for the Markov process, the n-point distribution 
can be written as 

Pn(xi,ui] . . .;x n ,u n ) = P{x n ,u n \x n - 1 ,u n -i) . . .P(x 2 ,u 2 \xi 1 ui)p 1 (x 1 ,ui) 
= K^P^ +1 -x % ,u M ) n> 

il= 2 piK) 

On the same line, a complete statistical description of shocks in Burgers solution is obtained through the n-shocks 
distribution densities which factorize to 

/ n n"=i P2(xi+i - Xi,Hi,rn;[H + i,T}i + i) 
pn{xi,Vi,'Oi,...,x n ,(i n ,r] n ) = -— T — , n>3. (99) 

lli=2 Pl{tM,Vi) 



The distribution of ordered sequences of next neighboring shocks is obtained from (|99|) by omitting the function H 
in p 2 , Eq. (f48|) . Here, factorization follows simply from the Markov property of Brownian motion and the fact that 
multiple constraints of the form ( |l4| ) decouple. From the point of view of the hierarchy of kinetic equations that 
governs the dynamics of shocks this factorization corresponds to an exact closure of this hierarchy or to an exact 
propagation of chaos. This will be discussed in ||. 
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As far as the time dependence is concerned, it can be reintroduced via the basic transition kernel ( |27| ) , which should 
be computed with s v (y) replaced by s u {y)/t. Owing to the invariance of the Brownian measure under the change 
if>(y) — > t l / s ip(y/t 2 / 3 ), one immediately finds that K u (fj,x, ryi, fi2, t) — t~ l l 3 K „> (ij/i , Wi , U^Vo) where the variables 
are rescaled according to ^ = /i^t -2 / 3 , rj[ — r\it~ l l 3 , and v' — vt~ 2 l 3 . From (|35|), (|47| ) and (|4(i| ) this implies the 
transformation laws of the functions J, /, and H 

J(u;t) = t~ 1/3 J{v'), I(n,r];t) = r 1 1 (ji' , T)') , and H(x, v x , u 2 ,t) = t- 2/3 H(x', u[, u' 2 ) (100) 

where x' = xt~ 2/>3 . This leads to the time dependent distributions 

p n (x 1 ,u 1 ; . . .;x n ,u n ;t) = t n/3 p n {x' 1 ,u' l ; . . .;x' n ,u' n ), (101) 

with v! i = Uit 1 / 3 , and 

p„(xi,/zi,7?i; . . .;x n ,^ n ,j] n ;t) = t~ 5n/3 p n (x' 1 , fi[, . . . ; x' n , fjf n> rf n ), (102) 



To obtain (101), we recall that the distributions p n {x\, u\\ . . . ; x n , u n ) were calculated from those of the coordi nate s 
of the contact points Xi — £j. At time t ^ 1, one has Xi — & = Ujf introducing a Jacobian i™ included in (101) 
when expressing the distributions as functions of the Burgers field amplitudes itj. From there, one recovers the 
well-known time dependent behavior of some moments of the distributions e.g., the energy dissipation per unit of 
length (u 2 {x,t)) ~ t~ 2 l 3 , the average number of shocks per unit of length ~ i~ 2 / 3 , the average strength of a shock 

(/i/t)~t-l/3. 
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APPENDIX A: 



We justify in this appendix the equation (pq ) which replaces the sum on zeros of the Airy function by an integral 
in the complex plane. 

To evaluate Ai(w) on the branch w + (p) = —p + iax^fp, we start from the formula Ai(— w) = e i7r / 3 Ai(u;e i7r / 3 ) + 
e -™/3Ai(we- l7T / 3 ) H giving^ 



Ai(w + (p)) = e™/ 3 Ai(-w+(p)e"/ 3 ) + e^' 3 A\{-w+{p)e^' 3 ) 
1 



( -i ^(- W + (p)) 3 / 2 J + exp (i 2 -{-w+{p)) 3 ' 2 



/47f(w+(p)) 1 /4 

where we have used the asymptotic behavior (|5^ ) of the Airy function Ai(w) for |«;| — > oo, arg w ^ ir. As p — > oo 

„3 



(Al) 



(- W +(p)) 3 / 2 = (p - i ax^) 3 / 2 = p 3 / 2 - i ^azp - \a 2 x 2 ^p - ? a 3 ^ + O 



Upon inserting (A2) in (Al) one sees that 



1 



M(w+(p)) 



< C p:x exp -axp - 



a 3 x 3 
24 



(A2) 



(A3) 



with Cp^ x growing at most algebraically with p and x. Using Ai(w;*) = Ai*(ui) one has the same estimate on the 
branch w~{p). By a similar calculation one has also that, for fixed £, Ai(£ + w ± (p)) / Ai(w ± (p) remains bounded as 
p — ► oo . 



6 The formula enables to obtain the asymptotic behavior of the Airy function Ai(u>) when argui approaches 7r as it is the case 
for w ± (p), p — > oo. 
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Consider now the finite parabolic contour closed by a circular arc TZe with 9 close to tt. On this circular arc for 
large radius TZ 



AU( + TZe w ) ( TZe i9 \ lfi ( 2_ m , fl ,, /2 2_ , fl ,, /2 \ 
" ) exp (--( C + He*)*" + -(^)3/ 2 j 

7?p^ \ 1/4 / \ 

' exp f-CV^e ifl / 2 ) = 0(1) (A4) 



Ai(fte if >) 

~ (c : - 'R< ' 

as TZ — ► oo and | < < 7r. Since fri + > 0, (r 2 — -j=j > the factors exp (awi (r% + an d 

exp ^au>2 (r 2 — "^?)) decay exponentially fast when w\ and w 2 are on the contour C x or on the circular arc. One 

concludes that the integrals on the circular arcs vanish as TZ — > oo so that the sums in ( |86| ) can indeed be replaced by 
the contour integrals. 



APPENDIX B: 

The integral B(x) 

p OO 

B{x)= / dCe C "Ai(C + wi)Ai(C + w 2 ) (Bl) 



is the Laplace transform for a negative argument —x of the product /(C) = Ai(C + u'i)Ai(C + i«2) of two Airy functions 
(omitting w\ and w% from the notation). First the asymptotic behavior of B(x) is determined by the Laplace method 

B(x) ~ -^e* (x) , x^oo (B2) 
2V 71 " 

where 

Hx) = -~-( Wl +w 2 )--\nx- K \ x ' ■ (B3) 

From the property of the Airy function (|2^), /(C) verifies the 4th order differential equation 

/""(C) - (4C + 2w x + 2w 2 )f'(0 - 6/'(C) + (wi - w 2 ) 2 f(0 = 0. (B4) 



From (B4), one finds that its Laplace transform for negative arguments satisfies 

B'(x) - h{x)B(x) = g{x) (B5) 

where we remark that 

h(x) = $'(x) (B6) 

and with 

= |/(o) - 1/'(0) + ^ [(/"(0) - 2( Wl + W2 )/(0)] 

- ^2 [./""(°) - 2 /( ) - 2(«* + «*)/'(0)] . (B7) 



Eq. ( p35[ ) can be solved, using also the value ( p2\) for a; — > oo, 

B(x) =B 1 (x) -B 2 (x) (B8) 

with 

1 f°° 
B 1 (x) = —=e*W, B2(x)=e*W dye^g(y). (B9) 

Notice that when evaluated at the zeros of the Airy functions w\ — —u>kn w 2 = —LUk 2 , s(y) reduces to 

/ \ I Ai(-Lj kl )Ai(-u> k2 ) 
g{y) \ Wl =- Uhl , W2 =-u h2 = ^ • ( mo > 
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APPENDIX C: 



We consider the multiple integral H\{x) ( p9[) and show first that it is absolutely convergent. On the contour 
w (p) = —p ± i ax^fp, < p < oo, we have 



W\ — W-2 

2a 2 ~Jx 



Pi - P2 



2a 2 v / i / 4a 2 



—^{ y /pl±^plf 



Hence, using (yfpi ± ^fpH) 2 < 2(p x + p 2 ) and (A3), the integrand ( P0[ ) is bounded by 
\h 1 (r,w 1 ,w 2 )\ < C Pu p 2iX exp 



12 



x exp < -ax(pi + p 2 ) - apivi + ap 2 v 2 - a 3 [ r + 



Pi - Pi V 
2a 2 Vi / ' 4 



<c 



: exp 



12 



exp < -apiri - ap 2 r 2 - a [r 



Pi ~ P2 
2a 2 y/x~ 



(CI) 



(C2) 



with C Pl}P2 ^ x increasing at most algebraically, showing that the integral ( p9[ ) converges absolutely 
To obtain the asymptotic behavior (B2h of Hi{x) we write the integration of hi(r, wi, w 2 ) over r as 



dr h 1 (r, wi,w 2 ) 



dr 



-%/xTx 



dr 



\fxri 



dr hi(r,wi,w 2 ). 



(C3) 



The first integration is readily performed to give J{vi)J(—v 2 ) (see Eq.(|9l|)) as J{v) ( |64| ) can be represented as an 
integral on any contour that encircles the zeros of the Airy function, in particular on C x . Thus it follows from ( |89| ) 
that 



Hxix)- J{V X )J{-V2) 



,5/2 



■ exp 



dr 



-a 3 v 3 + a 3 v\ 



dr 



*/xr 2 



{2mf 



dwi / dw 2 hi(rj,wi,w 2 ) 



(C4) 



Consider the contribution to (C4) where r > \fxr 2 and the branches of C x are w + (pi), w + {p 2 ). With (C2) this 
contribution is majorized by 



dw + (pi) 



dr / dwi / dw 2 h(r,wi,w 2 ) 



< exp 



12 



dpi 



dw+(p 2 ) 



dp 2 



C PuP2 , x exp < ~apin - ap 2 r 2 - a [r + 



dr / dpi I dp 2 

JO 



Pi - P2 



2a 2 y/x 



(C5) 



We split the p 2 integral into the domains < p 2 < a 2 yfxr and a 2 ^/xr < p 2 < oo. When < p 2 < a 2 yfxr, pi > 0, 
r > \fxr 2 = *Jx(xj2 — v 2 ) one has 



Pi ~ P2 V 

2a 2 y/x J ~ \2 ' 2a 2 y/x 



> I - 1 * 



> - > - 



> 



8 8 - 8 



(C6) 



where the last inequality holds for x large enough with c > 0, and thus 



exp < —a r 



Pi ~ P2 
2a 2 sfx 



< exp < —a 



ca 3 x 3 



(C7) 



On the other hand, when p 2 > a 2 yfxr : r > \/xr 2 , 



p 2 r 2 a 2 ryfxr 2 p 2 r 2 a 2 r\ p 2 r 2 3 
P2^2 > — 1 z > — z 1 -z— > — ^ h cx 



(C8) 



19 



where the last inequality holds for x large enough with c > 0. This leads to 

exp (— ap2f2) < exp (- 



ap2r 2 3 
acx 



(C9) 



The bounds (|C7|) and (C9) are introduced in (|C5|), the remaining r, pi, pi integrals are convergent and bounded with 



respect to x (except for a polynomial growth due to C PliP2jK and the line elements | dw ^ |= yl + The other 



contributions to (p4J) are treated in the same way. This leads to the result (92] 



APPENDIX D: 

We determine here the asymptotic behavior of G{x) (|9^) for large x. Starting from 



G{x)= J dz I dr^2 E &i r ' x ) 

Jo J-y/^ri kl >lk 2 >l 



(Dl) 



with 



Gk x ,k 2 {z,r;x) = 



exp 



V 12x a 



e a r exp -(wfej - WkJ' 



Aa{zx + ax 4 ) 



x exp 



r z \ I r : 
-auk! I n + —j= + - — - ) - aujk 2 I r 2 — + 



we define 



2ax 2 



F{x) = e aWiri+a " ir2 G(a;) 



2ax 2 



(D2) 



(D3) 



where r*i = x/2 + i/j, r 2 = x/2 — v 2 and — u>\ is the largest zero of the Airy function. We then decompose F(x) — 
F a (x) + Fb(x) + F c (x) according to the following splitting of the r integration range and the k±, fc 2 summations (for 
x large): 



oo pr 2 

dz / dre auJiri+aUir2 



F a (x) 
F b (x) 

F c (x) 

By dominated convergence, we immediately have that 



Gi,i(z,r;x) 



^ dz r d re a ^ +a ^ I E Gk u k 2 (z,r;x)-g hl (z,r;x)\ 

J ° J ~ r i \ki>lka>l J 



dz 



f\fxri 




\ dr + 




'1*2 





fci>l k 2 >\ 



/•oo / q2 g\ /*oo ^ 

lim F a (x) = I dz exp ( I / dr exp (— a 3 r 2 ) = — 

Jo V 4 / J-oo " a 



4^ 



4 J I ' <r - 

We then show below that F{,(x) and F c (x) vanish as x — > oo leading to the asymptotic behavior 

G{x) ~ l^L- Q <^+^\ (x - oo) 
a'' 



(D4) 
(D5) 

(D6) 

(D7) 
(D8) 



and thus to the behavior of H2(x), Eq . (B7[) . 

Since — r i < r < r 2 in the integral (|D5[) 7onc can choose x large enough so that r\ + -7= > rj. (1 — e) , r 2 — 7= > r 2 ( 1 — e) 
e > 0. Hence the fci, fc 2 term of the integrand in ( |D5| ) is less than 



(D9) 
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showing that the joint z, r in tegr als and fci, fc 2 summations converge. Moreover, since the term (fci, fc 2 ) = (1, 1) is 
absent from the integrand in (D5), there is at least one of the indices strictly greater than one. If both the indices 
are strictly greater than one, we can conclude that < Fb(x) < Cexp(— amin(ri,r 2 )(ti; 2 (l — e) — tends to zero 
exponentially fast as x — > oo provided that e < (cj 2 — lj\)/uj2 with — lj 2 the second zero of the Airy function. If one of 
the indices is equal to one, say fci = \,h%> 1, we have < F\,(x) < Cexp(— ar 2 (ot> 2 (l — e) — uji) + aer{) which tends 
exponentially to zero as x — > oo provided that e < (ll>2 — cji)/(1 + w 2 ). 

Consider now the integral in ( |D6| ) with r 2 < r < \fxr2- Since the factor exp (uj^ — ^fc 2 ) 2 4 a (za:+az 4 ) ) * s smaner 
than one, the ki, fc 2 summations are bounded by a product of X functions (|53|). Hence, for r > r 2 



< F c (x) < exp 



xe 



a 3 r^ 



dv 



dr exp 



lax 1 



r-2 



yfx 2ax 2 



For J m 



2ax* 



we use the bound X m + + 2^-^< Cexp ^-awi ^ri 



(D10) 



< Cexp(— a<jj\T\) 



since the argument becomes large as x — * oo, whereas for X ^r 2 — + 2a z a;2 j we use the bound (see the discussion 

-3/2 

since the argument can become small when r 



leading to Eq.@) 1 [r 2 



r)< C(r 2 



approaches the upper integration limit \fxr2- Thus 



< F c (x) < C 2 exp 



,3^2 



= C 2 4ax 5/2 exp 



,3^2 



aLO\T2 



dr exp 



r 2 exp 
dr 1 



dz 



) 2 



r2 



2ax* 



3/2 



(Dll) 



The second line has been obtain ed by performing the z-integral and changing the integration variable r to r' = 
T2 — r/yfx. This last integral in (Dll) is finite uniformly with respect to x so that with r 2 = x/2 — V2 the bound 
( Dll| ) tends to zero in a Gauss ian way as x — > oo. These last arguments can be reproduced to show that the integral 
with —y/xri < r < — n in Eq.(D6) tends to zero. 
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